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PROBLEM  STATEMENT  AND  MOTIVATION 

The  focus  of  this  research  effort  is  the  study  of  structural  dynamics  with  parameter  uncer¬ 
tainties.  The  motivation  for  this  study  rests  with  the  need  to  understand  the  dynamics  and 
control  of  large  space  structures  such  as  the  Space  Station.  The  parameter  uncertainties 
include  structural  mass,  damping,  and  stiffness,  and  environmental  loading. 

A  basis  for  this  study  comes  from  the  realization  that  numerous  uncertainties  exist  in  the 
modeling  and  analysis  of  large  space  structures.  These  include; 

1.  deployment  of  structural  components;  since  it  will  be  impossible  to  bring  into  orbit 
a  complete  space  structure,  components  will  be  carried  into  orbit  on  the  Space  Shuttle, 
then  deployed  and  constructed  in  a  manner  that  depends  on  the  specific  component.  De¬ 
ployment  dynamics  in  general  will  involve  uncertain  initial  conditions  and  may  depend  on 
environmental  factors  such  as  solar  radiation.  If  the  structure  to  be  deployed  is  of  the 
truss  type,  then  deployment  dynamics  (trsmsient)  must  account  for  uncert^linties  in  nodal 
coordinates  and  the  properties  (such  as  damping)  of  the  joints. 

2.  loading  uncertainties;  loads  on  the  space  station  will  be  due  to  either  unpredicted 
sources  or  to  known  but  very  complex  sources.  These  loads  will  be  transient  and  best 
modeled  as  random  processes.  Examples  of  such  loads  are  thermal  stresses,  dynamics  of 
deployment  and  operation,  and  impulses  due  to  impacts. 


3.  material  properties;  while  it  is  expected  that  the  properties  of  materi2ils  to  be  used  in  ~ 
space  station  design  will  be  well  known,  certain  structu^^ll  types,  such  as  composites,  may 
be  best  modeled  in  a  probabilistic  fashion.  Furthermore,  it  is  very  important  to  consider 
and  predict  fatigue  and  fracture  lives  of  Space  Station  components.  These  are  inherently 
probabilistic  problems. 

4.  extreme  behavior;  to  understand  the  limits  of  Space  Station  behavior  we  must  determine 

quantities  such  as  extreme  value  statistics  and  thus  able  to  establish  confidence  bounds  on  ,/ 
undesirable  behavior.  a..  l.  m  / 
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5.  parameter  estimation  and  system  sensitivity;  since  much  of  our  understanding  of  struc¬ 
tural  mechanics  is  derived  from  experimental  studies,  it  is  imperative  that  confidence 
bounds  be  placed  on  parameter  values  and  that  system  sensitivity  to  parameter  deviations 
be  established. 


TECHNICAL  APPROACH 

The  engineering  commimity  has  only  recently  become  convinced  of  the  necessity  to  account 
for  imcertainties  in  aii2ilysis  and  design.  The  use  of  probabilistic  modeling  techniques 
and  statistical  data  analysis  has  become  an  important  area  for  engineering  Research  and 
Development.  However,  to  date,  only  relatively  simple  structures  can  be  analyzed  where 
structural  pareuneters  and/or  environmental  loads  have  significant  random  components. 

Various  techniques  have  been  attempted  to  analyze  structures  where  significant  uncertain¬ 
ties  have  been  associated  with  mass,  damping,  or  stiffness.  These  techniques  were  based  on 
major  simplifying  assumptions.  In  the  current  effort,  an  “iterative/  decomposition”  tech¬ 
nique  was  adopted  that  permits  the  analyst  to  proceed  without  mzdcing  initial  assumptions 
or  simphfications.  In  this  manner,  a  most  general  approach  can  be  attempted  and  insur¬ 
mountable  mathematical  difficulties  can  be  addressed  as  they  appear  in  the  analysis.  Thus, 
needless  simplifications  may  be  avoided.  The  formal  technique  is  not  new,  but  its  previous 
applications  have  been  implemented  with  a  priori  simplifications  and  without  numerical 
results  for  verification.  The  present  effort  has  been  worked  through  to  numerical  results 
which,  by  way  of  a  Monte  Carlo  simulation,  will  be  verified  within  prescribed  error  bounds. 


SUMMARY  OF  RESULTS 

We  have  modeled  a  simple  structure  with  random  process  stiffness  and  environmental 
loads.  Two  physical  situations  which  may  be  modeled  in  such  a  feishion  are  structures  that 
have  parametric  loads  (pin-ended  column  with  axial  random  excitation),  and  situations 
where  the  structural  material  may  decay  in  a  random  meinner. 

Two  cases  have  so  far  been  considered  with  the  iterative/decomposition  technique: 

i)  non-white  (colored)  noise  random  processes,  and 

ii)  white-noise  random  processes. 

The  first  is  physically  more  realistic  and  much  more  difficult  mathematicjilly.  This  most 
general  case  was  undertaken  first  and  initial  analytic/numerical  (using  symbolic  manip¬ 
ulation  code  MACSYMA)  results  indicate  that  stiffness  parameter  uncertainties  have  a 
significant  effect  on  the  structural  dynamics.  Recently,  by  designing  a  Monte  Carlo  simu¬ 
lation  model  of  the  structure  and  environment,  we  have  been  able  to  estimate  errors  due 
to  approximations  in  (i).  These  compau'isons  with  simulation  have  permitted  us  to  sim¬ 
plify  the  approximate  technique  by  limiting  the  expansion  to  a  one-term  “correction”  to 
account  for  system  uncertsdnties. 
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The  second  case,  with  the  assumption  of  white-noise  random  processes,  is  much  simpler 
analytically,  and,  in  fact,  it  is  possible  to  obtain  exact,  closed-form  solutions. 

The  implications  from  these  results  are  the  following: 

(i)  parameter  \mcertainties  can  significantly  alter  structural  response  smd,  equally 
imp>ortant,  can  result  in  unstable  structural  behavior;  bounds  for  stable  behavior  have  to 
be  established  where  parametric  excitation  exists; 

(ii)  ranges  for  parameter  values  can  be  established  to  assure  not  only  stability,  but 
also  error  bounds. 

Additional  studies  were  performed  regarding  stochastic  stability;  this  has  been  accom¬ 
plished  in  addition  to  contracted  tasks  because  it  had  been  determined  that  such  work  weis 
needed  to  provide  a  more  complete  dimension  to  the  research  effort. 

Since  algebreiic/computational  codes  such  as  MACSYMA  are  relatively  new  to  such  work, 
it  was  decided  to  gather  our  experience  and  collaborate  with  two  academics  (DiMaggio 
and  Elishakoff)  to  issue  a  paper  summarizing  our  work. 


THIRD  YEAR  WORK 

The  Monte  Carlo  results  make  it  possible  to  consider  extensions  of  the  model  to  more 
complicated,  but  specific,  types  of  structures.  In  this  manner,  the  analysis  becomes  less 
general  when  specific  types  of  structures  with  more  narrowly  defined  uncertainties  are 
studied. 

We  have  concentrated  our  third  year  efforts  on  problems  of  stochaistic  stability,  stationarity 
of  stochastic  systems,  extending  our  understanding  of  this  expansion  method  to  leu’ge 
systems  by  way  of  “stochastic  finite  element  methodologies” ,  and  applying  the  method  to 
structures. 

Appendices  are  included  which  demonstrate  the  importance  of  incorporating  structural 
uncertainties  in  such  analyses.  Appendix  A  is  a  group  of  figures  where  only  the  input  is  a 
rjuidom  process.  These  curves  demonstrate  the  correctness  and  effectiveness  of  the  Monte 
Carlo  simulations  that  have  been  performed  for  purposes  of  model  verification.  Appendix 
B  is  a  group  of  curves  where  the  system  stiffness  parameter  K(<)  is  2issumed  to  be  positive- 
definite.  Here,  we  check  the  limits  of  the  derived  approximate  method  as  the  variance  of 
the  stochastic  parameter  is  increased  to  very  lau-ge  values  (1).  Appendix  C  is  a  group 
of  figures  that  demonstrate  the  effects  of  system  {K{t))  randomness  on  system  response. 
Two  facts  to  note  are  that  (i)  there  can  be  a  major  difference  in  response  variance,  and  (ii) 
there  can  be  major  differences  in  peak  “energy”,  i.e.,  the  peak  Spectral  Density  value  at 
system  natural  frequency  =  1))  will  differ  by  orders  of  magnitude.  All  these  issues 

have  been  discussed  in  detail  in  our  papers  and  reports. 

On  the  next  page  is  a  list,  by  year,  of  our  research  reports,  presentations,  and  published 
refereed  papers. 
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RANDOM  VS.  DETERMINISTIC  SYSTEM. 


VARIANCE  OF  1C  -  0.50,  BETA  -10.00,  ALPHA  -  1.00 
VARIANCE.  DETERMINISTIC  -  2.73,  RANDOM  -***•* 


—I - , - , - 1 - , — 

10.0  20.0  30.0  400  60.0 

DIMENSIONLESS  TIME. 


RANDOM  VS.  DETERMINISTIC  SYSTEM. 


VARIANCE  OF  K  -  0.40.  BETA  -10.00,  ALPHA  -  1.00 
PEAK  AT  1.0,  DETERMINISTIC  -25.00  .  RANDOM  -***•« 


RANDOM  K 


DETERMINISTIC  K 


AUTOCORRELATION. 

20.0  *10.0  0.0  10.0  20.0  30.0  40.0 


RANDOM  VS.  DETERMINISTIC  SYSTEM. 


VARIANCE  OF  K  -  0.40.  BETA  -=10.00,  ALPHA  -  1.00 
VARIANCE.  DETERMINISTIC  -  2.73.  RANDOM  -34.43 


RANDOM  VS.  DETERMINISTIC  SYSTEM. 
VARIANCE  OF  K  -  0.30.  BETA  -10.00.  ALPHA  -  1.00 
PEAK  AT  1.0.  DETERMINISTIC  -25.00  .  RANDOM  -***•< 


RANDOM  K 


DETERMINISTIC  K 


AUTOCORRELATION. 


RANDOM  VS.  DETERMINISTIC  SYSTEM. 
VARIANCE  OF  K  -  0.30,  BETA  -10.00.  ALPHA  -  1.00 
VARIANCE.  DETERMINISTIC  -  2.73,  RANDOM  -  9.98 


o 


RANDOM  VS.  DETERMINISTIC  SYSTEM. 


VARIANCE  OF  K  -  0.30.  BETA  -10.00,  ALPHA  -  0.10 
PEAK  AT  1.0,  DETERMINISTIC  -  4.95  .  RANDOM  -58.35 


AUTOCORRELATION. 


RANDOM  VS.  DETERMINISTIC  SYSTEM. 


VARIANCE  OF  K  -  0.30.  BETA  -10.00,  ALPHA  -  0.10 
VARIANCE.  DETERMINISTIC  -  1.46.  RANDOM  -  5.99 


RANDOM  K 


DETERMINISTIC  1C 


A 


SPECTRAL  DENSITY. 

BO.O  tOO.O  1B0.0  200.0 


RANDOM  VS.  DETERMINISTIC  SYSTEM. 


I 


VARIANCE  OF  K  -  0.36.  BETA  -10.00.  ALPHA  -  0.10 
PEAK  AT  1.0.  DETERMINISTIC  -  4.95  ,  RANDOM 


RANDOM  K 


DETERMINISTIC  K 


j 


RANDOM  VS.  DETERMINISTIC  SYSTEM. 
VARIANCE  OF  K  -  0.35,  BETA  -10.00,  ALPHA  -  0.10 
VARIANCE.  DETERMINISTIC  -  1.46,  RANDOM  -  8.81 


RANDOM  K 


DETERMINISTIC  K 


SPECTRAL  DENSITY. 

60.0  100.0  160.0  200.0  260.0  300.0 


RANDOM  VS.  DETERMINISTIC  SYSTEM. 
VARIANCE  OF  K  -  0.40.  B^A  -  1.00.  ALPHA  -10.00 
PEAK  AT  1.0,  DETERMINISTIC  -25.00  .  RANDOM  -95.65 


RANDOM  K 


DETERMINISTIC  K 


\\  h  IIM.INCKH  \S>(h;|  M1>> 
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RANDOM  VS.  DETERMINISTIC  SYSTEM. 


VARIANCE  OF  K  -  0.20.  BETA  -  1.00,  ALPHA  -10.00 
VARIANCE.  DETERMINISTIC  -  2.50.  RANDOM  -  3.86 


RANDOM  K 


DETERMINISTIC  K 


RANDOM  VS.  DETERMINISTIC  SYSTEM. 


VARIANCE  OF  K  -  0.40.  BETA  -  1.00,  ALPHA  -  1.00 
PEAK  AT  1.0,  DETERMINISTIC  -25.00  ,  RANDOM  -*•** 


o 

6 


w  kii)I.in(;kh  vssoci  \tks 


vwjfwy, 


AUTOCORRELATION. 

0.0  1.0  2.0  3.0  4.0 


W  F  IDl  INf.KH  \s>(l(  1M  K> 


I 

•:A^ 


J 


AUTOCORRELATION. 


!*•••*  ill 


k!i*^i*j 


SPECTRAL  DENSITY. 

p.O  S.O  10.0  16.0  20.0  26.0  30.0  36.0  40.0 


RANDOM  vs.  DETERMINISTIC  SYSTEM. 


VARIANCE  OF  1C  -  0.30.  BETA  -  1.00.  ALPHA  -  0.10 
PEAK  AT  1.0.  DETERMINISTIC  -  4.95  .  RANDOM  -31.43 


AUTOCORRELATION 


RANDOM  VS.  DETERMINISTIC  SYSTEM. 


VARIANCE  OF  K  -  0.30.  BETA  -  1.00,  ALPHA  -  0.10 
VARIANCE.  DETERMINISTIC  -  1.46,  RANDOM  -  7.09 


W  K1I)1.IN(.KH  \SS(»(.l  \TKS 


SPECTRAL  DENSITY.  •lOf 

p.O  6.0  10.0  16.0  20.0  26.0  30.0 


RANDOM  vs.  DETERMINISTIC  SYSTEM. 
VARIANCE  OF  K  -  1.00.  BETA  -  1.00,  ALPHA  «  0.10 
PEAK  AT  1.0,  DETERMINISTIC  -  4.95  .  RANDOM 


AUTOCORRELATION. 

-2.0  -1.0  0.0  1.0 


RANDOM  VS.  DETERMINISTIC  SYSTEM. 
VARIANCE  OF  K  -  1.00.  BETA  -  1.00.  ALPHA  -  0.10 
VARIANCE.  DETERMINISTIC  -  1.46.  RANDOM 


"to  9 


W  (  Ilt|  IN(.KK  \»<  M  I  Ml  ^ 


AUTOCORRELATION. 


AUTOCORRELATION. 


K 


W  Kill!  I\(.KK  \SNi(  I  \TI> 


RANDOM  VS.  DETERMINISTIC  SYSTEM. 


VARIANCE  OF  K  -  0.15.  BETA  «  0.10,  ALPHA  «  1.00 
PEAK  AT  1.0,  DETERMINISTIC  -25.00  ,  RANDOM  -23.01 


1  RANDOM  K 


DETERMINISTIC  K 


— I - , - 1 — 

3.0  4.0  6.0 

ANGULAR  FREQUENCY. 


— 

e  0 


-I - 1 

7.0  8.0 


RANDOM  VS.  DETERMINISTIC  SYSTEM. 
VARIANCE  OF  1C  »  0.15.  BETA  *  0.10,  ALPHA  -  1.00 
VARIANCE.  DETERMINISTIC  -  2.73.  RANDOM  -  8.80 


RANDOM  K 
DETERMINISTIC  K 


1 

1 

3 

i 

W  KII)|.IN(.KH  \SS()(;i  \  Y\ 

RANDOM  VS.  DETERMINISTIC  SYSTEM. 


VARIANCE  OF  K  -  0.10.  BETA  «  0.10,  ALPHA  «  1.00 
PEAK  AT  1.0,  DETERMINISTIC  -25.00  .  RANDOM  -14.62 


AUTOCORRELATION. 

900.0  0.0  600.0  1000.0  1600.0  2000.0 


RANDOM  VS.  DETERMINISTIC  SYSTEM. 


VARIANCE  OF  K  -  0.20.  BETA  -  0.10,  ALPHA  -  1.00 
VARIANCE.  DETERMINISTIC  -  2.73.  RANDOM  -***•* 


o 

6 

o 

o 


SPECTRAL  DENSITY. 

1000.0  2000.0  3000.0  4000.0  SOOO.O  6000.0  7000.0  8000.0  9000.0 


RANDOM  VS.  DETERMINISTIC  SYSTEM. 


VARIANCE  OF  K  -  0.20,  BETA  -  0.10,  ALPHA  -  1.00 
PEAK  AT  1.0,  DETERMINISTIC  -25.00  ,  RANDOM  -****• 


RANDOM  VS  DETERMINISTIC  SYSTEM. 


VARIANCE  OF  K  -  0.15  BETA  -  0.10  ALPHA  -  0.10 
VARIANCES.DETERMINISnC  -  1.46  RANDOM  -15.81 


RANDOM  K 


DETERMINISTIC  K 


10.0 


20.0  30.0  40.0 

DIMENSIONLESS  TIME. 


60.0 


eo.o 


SPECTRAL  DENSITY. 

0.0  10.0  20.0  300  40.0  60.0  60.0  70.0  60.0  60.0 


w  kii>i.in(;kk  associ  atks 


1 


I -i'S? 


sv 


SPECTRAL  DENSITY. 

10.0  18.0  20.0  28.0  30.0 


W  KIDI.IM.KH  \SSO(.l  VT  KS 


•7>S5 


Gsa- 


RANDOM  VS.  DETERMINISTIC  SYSTEM. 


VARIANCE  OF  K  -  0.10,  BETA  -  0.10,  ALPHA  -  0.10 
VARIANCE.  DETERMINISTIC  -  1.46,  RANDOM  -  2.57 


RANDOM  K 


DETERMINISTIC  K 


“I — 
10.0 


I  1  I 

20.0  30.0  40.0 


DIMENSIONLESS  TIME. 


60.0 


ei 


